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I. INTRODUCTION 

Low dimensional systems have been the subject of in- 
creasing interest over the past decades, traditionally in 
condensed-matter physics, and recently in the framework 
of quantum gases. There, recent advances in atom waveg- 
uide technology Q, |^ Q |^ , the realization of quantum 
gases in very anisotropic traps and loading Bose- 
Einstein condensates (BEC) on optical lattices 
allowed experimentalists to obtain a rich variety of sys- 
tems where the reduced dimensionality rules the physics. 
Particularly interesting is the case in which the quantum 
dynamics of the system becomes quasi-one-dimensional. 
In that case it was shown ^J, 0, that in certain 
regimes of large scattering length, low densities, and low 
temperatures, bosons behave as impenetrable particles 
known as hard-core bosons (HCB's). 

The one-dimensional (ID) homogeneous gas of HCB's 
was introduced by Girardeau |l4| , and a one-to-one corre- 
spondence between ID HCB's and spinless fermions was 
established. It was shown later by Lenard and by Vaidya 
and Tracy 0, 0| that in homogeneous space, the ID 
HCB gas does not exhibit true condensation at zero tem- 
perature since there is only off-diagonal quasi-long-range 
order (one-particle correlations decay as a power law), 
leading to an occupation of the lowest effective single- 
particle state (the highest occupied one) \/Nb, where 
Nh is the total number of HCB's. 

Recently, theoretical studies have been focused on the 
ground-state properties of ID HCB gases in the presence 
of harmonic potentials jHlil El iH ISl IH , which are 
required in the experiments with quantum gases to keep 
the particles confined. Using the Fermi-Bose mapping, 
quantities like density profiles, momentum profiles, the 
natural orbitals (effective single-particle states), and the 
one-particle density matrix have been calculated for a 
finite number of particles. Some results were extended 
to the thermodynamic limit {N^, oo) 0,1121 . It has 
been found that similarly to the homogeneous case the 
occupation of the lowest natural orbital is ~ ^/Nl [20ll2ll |. 
so that only a quasicondensate develops in the trap. 

Another case of interest is given by HCB's on a lat- 
tice. These systems have been realized experimentally 



very recently [l^. In the periodic case, the ID HCB 
Hamiltonian can be ma ppe d into the ID XY model of 
Lieb, Schulz, and Mattis [2J|. This model has been exten- 
sively studied in the literature. At zero temperature the 
asymptotic behavior of the correlation functions is known 
j2J, |2y, 1231 . In the presence of a confinin g po tential, the 
systems realized experimentally in Ref. |23| . the situa- 
tion is more complicated. This is because the trapping 
potential is equivalent to the addition of a space-varying 
transverse field to the XY Hamiltonian, for which in gen- 
eral analytical results are not available. 

In this work we study ground-state properties of ID 
trapped HCB's in a lattice. The lattice opens new possi- 
bilities for engineering states that cannot be obtained in 
continuous systems. In particular, for HCB's it is possi- 
ble to create pure Fock states when the occupation in 
some regions of the system reaches the value n = 1, 
such that coherence is lost there, and these sites decouple 
from the rest of the system. In addition, the properties 
of the systems without the lattice can be recovered in 
the low-density limit in a lattice. We follow a numerical 
approach, based on the Jordan- Wigner transformation 
|28j . that allows us to calculate exactly the one-particle 
Green's function and analyze arbitrary confining poten- 
tials. This approach was recently introduced in Ref. (2^ 
to study ground-state properties, and generalized in Ref. 
|30l | for the study of the nonequilibrium dynamics. 

We concentrate here on the behavior of the off-diagonal 
part of the one-particle density matrix, the momentum 
distribution function, and the natural orbitals. These 
quantities have a different behavior for HCB's when com- 
pared with spinless fermions to which the HCB Hamil- 
tonian is mapped by means of the Jordan- Wigner trans- 
formation. We find that the one-particle density matrix 
{pij) decays as a power law \xi — Xjj"^/^ for large dis- 
tances, like in the periodic case jl^, irrespective of the 
confining potential chosen. This is valid even when por- 
tions of the system reach occupation ni = 1, such that 
coherence is lost there. The power law above is shown 
to determine the scaling of the occupation of the lowest 
natural orbital in the thermodynamic limit. This scal- 
ing and its finite-size corrections are also studied for ar- 
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bitrary powers of the confining potential. In addition, 
the low-density limit in the lattice is analyzed in detail 
and the results obtained are compared with the ones pre- 
sented in Ref. for continuous systems. 

The exposition is organized as follows. In Sec. II we 
describe the numerical approach used. In Sec. Ill we 
discuss the properties of HCB's on periodic lattices. Sys- 
tems confined in harmonic traps are analyzed in Sec. IV. 
In Sec. V we generalize to arbitrary powers of the confin- 
ing potential the results obtained for the harmonic case. 
Finally, the conclusions are given in Sec. VI. 



II. EXACT APPROACH 

In the present section we describe in detail the exact 
approach we follow to study ID HCB's in the presence of 
a lattice I23. The HCB Hamiltonian can be written as 



H 



-t 



h.c. 



y^^x^n,, (1) 



where 1^'^^^) HCB ground-state wave func- 

tion, and 1^'^) is the equivalent noninteracting fermion 
ground-state wave function. In addition, we denote 



1*^) = 




(6) 



The ground-state wave function of the equivalent 
Fermionic system can be obtained by diagonalizing Eq. 
(IH) (the properties of such systems were analyzed in Ref. 
|3l|), and can be written in the form 



n 



N 



P^sfl |0), 



(7) 



<T=1 



with the addition of the on-site constraints 



^t2 



(2) 



These constraints on the creation (&|) and annihilation 
(bi) operators avoid double or higher occupancy. Notice 
that the brackets in Eq. Q apply only to on-site anticom- 
mutation relations, while = for i 7^ j, i.e., they 

commute on different sites |24| . In Eq. the hopping 
parameter is denoted by t, and the last term describes an 
arbitrary confining potential, with power a and strength 
Va', fii = blbi is the particle number operator. 

In order to obtain the ground-state properties of 
HCB's, we use the Jordan- Wigner transformation [28|. 



(3) 



13=1 13=1 

which maps the HCB Hamiltonian into the noninteract- 
ing spinless fermions Hamiltonian 



h.c. 



(4) 



where /| and fi are the creation and annihilation opera- 
tors for spinless fermions and nj = f- fi is their particle 
number operator. This means that HCB's and fermions 
have exactly the same spectrum. The nontrivial differ- 
ences between both systems appear in the off-diagonal 
correlation functions as shown below. 

Using the mapping above, the one-particle Green's 
function for the HCB's can be written in the form 



-1 



HCB/ 



i-1 



/3=1 
(*Fl*f), 



7=1 



with Nf the number of fermions (Nf = Nf,), N the num- 
ber of lattice sites, and the matrix of the components P 
is given by the lowest Nf eigenfunctions of the Hamilto- 
nian, 



P21 



P12 
P22 



\Pni Pi 



N2 



PlNf \ 
P2N, 



PNNf I 



(8) 



In order to calculate |^;^) (and )) we notice that 



n = n [1 - 2 !iu 



13=1 



0=1 



(9) 



Then, the action of Yl 



-i^ffffii Fermionic 
ground state [Eq. ITJ] generates only a change of sign 
on the elements P^s for cr < i — 1, and the further cre- 
ation of a particle at site i implies the addition of one 

1 and all the 

B\ 



column to P with the clement R 



Nf + l 



others equal to zero (the same can be done for IvE*]?)) 
Then 1*^) and |*f ) can be written as 



N' 



f N 



5=1 a=l 



N', 



f N 



1*^) = n E 10)' 



(10) 



5=1 a=l 



where P' and P' are obtained from P changing the 
(5) required signs and adding the new column (iV'^=A^/-|-l). 
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The Green's function Ga is then calculated as 



N 



AT', 



N 



5=1 5-=l 



5=] 
N 

E 



(5=1 <T=1 



p/A p/, 

^(Ti 1 • • • 



.lA 

a„,N> 



x(0|/., 



det 



f V- 

t , 



/IjO) 



(11) 



where the following identity was used 



Ai-A 



(12) 

with e ^ "j" the Levi-Civita symbol in N'^ dimensions, 
the indices A have values between one and N'^. 

The Green's function Hll|l is then evaluated numeri- 
cally, and the one-particle density matrix is determined 
by the expression 



P^^, = (^I^j) = G^J + S,j (1 - 26,, 



(13) 



An alternative approach to the one presented here to 
calculate the one-particle density matrix of HCB's con- 
fined on optical lattices, using the Jordan- Wig ner trans- 
formation, was followed by Paredes et al. |23|| . In their 
work the elements of the one-particle density matrix were 
evaluated as Toplitz determinants of matrices with sizes 
up to [N — 1) X {N — 1). Within our approach only de- 
terminants of (A^6 + 1) X [Nb + 1) matrices are evaluated. 



III. HARD-CORE BOSONS ON PERIODIC 
SYSTEMS 

In this section we analyze the properties of HCB's in 
ID periodic lattices. In this case the Hamiltonian can 
be written as 



H ^-t 



+1 



h.c 



(14) 



This Hamiltonian, with the additional constraints Q, is 
particle-hole symmetric under the transformation hi = 
bl, hi = bi. The operators /ij and hj are the creation and 
annihilation operator for holes. The previous symmetry 
has important consequences since it implies that the off- 
diagonal elements of the one-particle density matrix for 
Nb HCB's [pij (Nb)] and for N - Nb HCB's [py {N ~ Nb)] 
are equal. Only the diagonal elements change, and they 
satisfy the relation pii{Nb) = 1 — pii{N — Nb). 

In the momentum distribution function (rjfe), the 
particle-hole symmetry leads to 

N-Nb' 



UkiNb) = n.kiN - Nb) + 1- 



N/2 



(15) 



so that an explicit dependence on the density appears. 
This behavior is different from the one of noninteract- 
ing spinless fermions, where the Hamiltonian is also 
particle-hole symmetric but under the transformation 
h[ = (—1)'//, h'^ ~ (— l)*/i. In this case the mo- 
mentum distribution function for fermions ['nj.{Nf)] and 
[nl(N — Nf)] are related by the expression 



iNf) = l-ni,^^{N^Nf), 



(16) 



where there is no explicit dependence on the density. 
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1: Momentum profiles for periodic systems with 1000 
lattice sites and occupations of 101 (a), 501 (b), and 799 (c) 
HCB's. 

In Fig. n we show the momentum profiles for systems 
with three different fillings in 1000 lattice sites. No- 
tice that we consider odd numbers of particles, which 
allows using periodic boundary conditions [23|. (In the 
case of an even number of particles antiperiodic bound- 
ary conditions are required •) The peak structure in 
the momentum distribution function (Fig. ^ refiects the 
Bosonic nature of the particles, and is in contrast with 
the structure of the momentum distribution function for 
the equivalent noninteracting fermions. On increasing 
the number of HCB's two effects can be seen: (i) up to 
half filling the value of rifc^g increases, and it starts to 
decrease when the number of particles exceeds half fill- 
ing; (ii) the population of high momenta states is always 



4 



increasing, showing that on increasing the density in the 
system, HCB's become more locahzed. The latter can be 
understood due to the impenetrability property of HCB's 
and the ID character of the system. 



A. Off-diagonal correlations and Tife^o 

The formation of the peak in the momentum distri- 
bution function of the HCB's is due to the off-diagonal 
quasi-long-range correlations present in the one-particle 
density matrix. As is shown in Fig. [51 they decay as 
Px = A{p)/ \/ xja (x = \xi — Xj\). Ai^p) is a function 
of the density (p = N^/N) in the system. A trivial cal- 
culation shows that in the thermodynamic limit, keeping 
the density constant, n^^Q ^ ^/-^J^ij scales 
as 

nk^o^B{p)^b^Cip)VN. (17) 

where B{p) and C{p) are functions of the density. The 
previous result is valid to a good approximation for fi- 
nite systems (for occupations higher than 100 HCB's), 
as shown in the inset in Fig. [21 where we plot nk=o vs 
Nb for different systems at half filling. The straight line 
displays a v^iV^ behavior. 
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FIG. 2: One-particle density matrix vs x/a for periodic 
systems with N = 1000 and p = 0.5 (dotted line), p = 0.211 
(dashed-dotted line), p = 0.071 (dashed line). Thin contin- 
uous lines correspond to power laws xja. The inset shows 
Tife^o vs A^f, for systems at half filling (O)) Ihe line exhibits 
\fNl behavior. 

A global picture of the occupation of the state with 
zero momentum as a function of the density is shown in 
Fig. [3 There we plot 7iJ,^o ~ ?^fc=o/v^ vs the density, 
and compare systems with 1000 and 300 lattice sites. The 
comparison shows that already for these system sizes the 
finite-size corrections are small. Actually, they are found 
to be given by 

nfc=o/Vlv = C(p)-i?(p)/ViV, (18) 

where £'(/o) can be positive or negative depending on the 
density. The inset in Fig.[21shows M^^^ = C{p) — n'f.^Q vs 



N~'^/'^ for half filled systems. The straight line displays 
the result of our fit for C{p = 0.5) and D{p = 0.5), and 
confirms the validity of Eq. (|18|) . The same functional 
form of the finite-size corrections was obtained for the 
continuous case |2ll |. which corresponds to systems in a 
lattice where the average interparticle distance is much 
larger than the lattice constant, i.e., in the limit p = 



1 I 1 1 1 1 1 1 1 r 




FIG. 3: Normalized occupation of = {n'^.^Q) as a function 
of the density in periodic systems with N = 1000 (thick con- 
tinuous line) and A'^ — 300 (dashed line) . The thin continuous 
line shows a fit of Eq. 12111 . with 6 — 1.51 and c — 1.18, to our 
numerical results for TV — 1000. The inset displays Sn'i.^Q vs 
TV"^/^ (see text) for half filled systems (O)) the straight line 
shows the result of our fit to Eq. (dHJ. 



B. Low-density limit in the lattice 

In what follows we study the corrections introduced 
by the lattice to the occupation of the lowest momentum 
states in the continuous case. In a lattice at low mo- 
menta, the spectrum is quadratic in k and Sk ^ 1/N, so 
that the ratio between the level spacing and the band 
width reduces proportionally to 1/A^^ when N is in- 
creased. We find that the low-momentum occupations 
approach their values in the continuous case in exactly 
the same way, 

A,(iVfc,iV)==A„(iV,)-:^^, (19) 

where A,, {Nb,N) = nk=2TT'n/Na (-^fe) is the occupation of 
a momentum state rj when there are Nb HCB's in a sys- 
tem with A'' lattice sites. A,, (Nb) is the occupation of the 
momentum state 77 in the continuous case, and Eri (Nb) 
for a given 77 is a function of Nb- 

We obtain A^ (Nb) and E^i (Nb) for the three lowest 
momenta states analyzing systems with fillings up to 501 
particles and sizes up to 2000 lattice sites. In Fig. ^ 
we compare our results for A^ (N h) w ith the ones in the 
continuous case presented in Ref. |23| [see Eqs. (56)-(58) 
there]. The agreement is excellent. In the inset we plot 
AA = A^ (A^h) — A,, {Nb, N) also for the three lowest mo- 
mentum states as a function of 1 /N'^ confirming Eq. H19|l . 
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FIG. 4: Extrapolated values of the first (0)i second (A), and 
third (O) momentum-state occupations in the continuous case 
(A,,). The lines following the data correspond to the results 
obtained in Ref. |2l| . The inset shows AA (see text) vs 
also for the first (0)i second (A), and third (O) momenta 
states in systems with 501 HCB's and sizes between 1000 and 
2000 lattice sites, the straight lines are the result of our fits. 

Surprisingly, we find that Eq. H19|l works extremely well 
up to very high densities (n ~ 0.8). This feature allows 
one to understand the dependence of the lowest natural 
orbital occupation on the density (Fig. O as follows. 

The dependence of Aq (A^b) and Eq (Ni,) on Ni, can be 
determined given Eq. ^ to be 

Ao (Nb) = b v/a^, Eo im) = c nI'\ (20) 

where b and c are two constants that can be obtained 
numerically. (Both relations above were confirmed by 
our numerical results.) Then Eq. H19|) may be rewritten 

as 

\Q{N,p)|^/N = ^p{b-cp'). (21) 

A fit of Eq. (|2I1) to our numerical results for N = 1000, 
with b = 1.51 and c = 1.18, is shown in Fig. O as a thin 
continuous line. As it can be seen, Eq. H21() describes 
very well the numerical results up to densities n 0.8, 
which were the densities up to which we had found that 
Eq. H19f) was a good approximation. 

C. nfe=o singularity in the thermodynamic limit 

To conclude this section on periodic systems, we ana- 
lyze in detail the low-momentum region of the momen- 
tum distribution function. The power -law decay of the 
one-particle density matrix {p^ ~ \j \fxja) implies that 
in the thermodynamic limit the momentum distribution 
hmction has a \k\~^ {P = 1/2) singularity at fc = [l6l |. 
For finite systems we find that, due to finite-size effects, 
the apparent exponent observed in the low-momentum 
region of Uk depends strongly on the number of particles 
in the system. In Fig. Owe show a log- log plot of the mo- 
mentum distribution of two systems at half filling with 
N = 100 and TV = 1000 lattice sites. A fit to a power 



law in the case of 100 lattice sites reveals an exponent 
P = 0.63. As the figure shows, the increase of the sys- 
tem size reduces the exponent of the power law only in 
the low-momentum region. We have studied systems at 
half filling fitting power laws from the lowest ten momen- 
tum states available (excluding fc = 0). The results are 
shown in the inset of Fig. [S] This inset shows that with 
increasing the system size the exponent (3 extrapolates 
very slowly to 0.5. Basically, we find that it extrapo- 
lates as a power law ~ jv-o.63_ rpj_^g latter exponent is 
a nonuniversal one since it depends on the way in which 
one fits power laws for the low- momentum region of n^. 
Up to the bigger system we calculated N = 1500, the 
exponent (3 reduced only up to 0.523. Our results imply 
that the experimental observation of the thermodynamic 
limit exponent (3 — 0.5 will be very difficult if not impos- 
sible considering finite-temperature effects |2t| . 




FIG. 5: (Color online) Log- log plot of the momentum distri- 
bution function of systems at half filling with A'^ — 100 (Q) 
and A' = 1000 (+) lattice sites. The straight line shows a 
power-law fit Uk ~ k~'^ (/3 = 0.63) to the lowest momentum 
states of the system with N — 100. The inset displays the de- 
pendence of the exponent (3 (fitted for half filled systems) on 
the system size. The straight line shows that it extrapolates 
to f3 = 0.5 for N ^oo. 



IV. HARD-CORE BOSONS CONFINED IN 
HARMONIC TRAPS 

We analyze in this section the case in which HCB's 
arc confined in harmonic traps. The Hamiltonian in this 
case is given by Eqs. ^ and Q with a = 2, and is not 
particle-hole symmetric like in the periodic case. In or- 
der to quantitatively characterize these systems we make 
use of the length scale set by the combination lattice- 
confining potential C = (^a/t) i a-nd the associated 
characteristic density p — N^/C,, both introduced in Refs. 
[31I Is^ . Is^ . Since the diagonal elements of the one- 
particle density matrix for HCB's and the equivalent 
fermions are equal [see Eq. Q for i = j], density prop- 
erties of both systems trivially coincide. We concentrate 
here in the quantities related to the ofF-diagonal corre- 
lations of the HCB's, like the momentum distribution 
function and the natural orbitals. 
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The natural orbitals {(f^) are defined as the eigenfunc- 
tions of the one-particle density matrix |34j . 



N 



(22) 



and can be understood as being effective single-particle 
states with occupations A,,. The lowest natural orbital 
(the highest occupied one) is considered to be the con- 
densate. In the periodic case the natural orbitals are 
plane waves, i.e, the eigenfunctions of Eq. H22|) are mo- 
mentum states (the lowest natural orbital is the state 
with momentum zero). On the contrary, since in the 
trapped case the translational invariance is broken, the 
natural orbitals occupations and the momentum distri- 
bution function do not coincide. 
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FIG. 6: Density (a)-(c) and normalized momentum distribu- 
tion function (d)-(f) for trapped systems with 1000 lattice 
sites, V2a^ = 3 x lQ~^t and occupations of 101 (a),(d), 401 
(b),(e), and 551 (c),(f) HCB's. 

In Fig. El we show three density profiles and their 
corresponding normalized momentum distribution func- 
tion for traps with 1000 lattice sites. Positions are nor- 
malized by the characteristic length and the nor- 
malized momentum distribution function is defined as 
"fc = {aK)Y.fj=ie-'''^"~'Kb\h,). Figure displays fea- 
tures similar to the periodic case. The normalized mo- 
mentum distribution function exhibits narrow peaks at 
fc = with nk=o initially increasing with the number of 
particles, and decreasing after certain filling when more 
particles are added to the system. Localization effects 
also start to be evident at high fillings, with an incre- 
ment of the population of high momenta states [Fig.|SJf )] . 



and the formation of Fock states (or a Mott insulating 
plateau) in the middle of the system [sites with 71^ = 1 
in Fig. EJc)] . Notice that in the latter case the forma- 
tion of the Mott insulator did not destroy the sharp peak 
observed at nk=a, which is associated to the superfluid 
phase surrounding the Mott plateau. 
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FIG. 7: Occupation of the natural orbitals for trapped sys- 
tems with 1000 lattice sites, V2a^ = 3 x W~^t, and occupa- 
tions of 101 (a), 401 (b), and 551 (c) HCB's. The insets show 
the occupation of the lowest eleven natural orbitals in each 
case. 

Results obtained for the occupations of the natural or- 
bitals, in the systems of Fig. El are presented in Fig. [7\ 
The occupations are plotted as a function of the orbital 
numbers ry, and they are ordered starting from the high- 
est one. The effects of increasing the filling in the system 
are similar to the ones observed in the momentum dis- 
tribution function. The occupation of the lowest natural 
orbital increases with increasing the number of particles 
up to a certain filling where it starts to decrease when 
more particles arc added. In contrast to the momentum 
distribution function, the natural orbital occupations ex- 
hibit clear signatures of the formation of the Mott insu- 
lating state in the middle of the system since a plateau 
with A = 1 appears [Fig. [Tfc)] , and the natural orbitals 
with A 7^ 1 become pairwise degenerated [inset in Fig. 
Uc)]. This can be easily understood since in the region 
where = 1, the one-particle density matrix is diagonal 
Pij ~ , so that it can be seen as a three block diagonal 
matrix with the block in the middle being an identity 
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matrix and the other two being identical. Thus the di- 
agonalization gives a group of eigenvalues = 1 and all 
the others pairwise degenerated. 
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FIG. 8: Profiles of the two lowest natural orbitals for trapped 
systems with: 1000 lattice sites, V2a? = Z x 10~^t, and oc- 
cupations of 101 (a),(b), 401 (c),(d), and 551 (e),(f) HCB's 
(continuous line); 300 lattice sites, V2a? = 3.3 x 10~*t, and 
occupations of 30 (a),(b), 121 (c),(d), and 167 (e),(f) HCB's 
(dashed line). 

Profiles of the two lowest natural orbitals, for the same 
parameters of Figs.|niand[71 are presented in Fig.|Sl Posi- 
tions in the trap arc also normalized by C, and the proper 
definition of the scaled natural orbitals is given by 



R 



1/2^') 



1/2 



(23) 



This scaling relation is valid only for the lowest natural 
orbitals and is meaningful for comparing systems with 
the same characteristic density. Notice that for a given 
p the occupied system size (L) is proportional to i.e, 
L = F{p)( with F{p) depending only on the value of p. 
In Fig.|Slwe have also plotted results for smaller systems 
(dashed line) but fulfilling the condition of equal charac- 
teristic density. The comparison between the plots shows 
that the scaling relation defined by Eq. (|23|l holds. (We 
find this to be valid to a good approximation for occupied 
system sizes larger than 100 lattice sites.) In Figs. |HIa) 
and (b) it can be seen that for low characteristic densi- 
ties the natural orbitals are similar to the ones of sys- 
tems without the optical lattice [20l I2H . On increasing 
the characteristic density the shape of the lowest natural 
orbital starts to change since its weight starts to reduce 



in the middle of the trap, as shown in Fig. |5^c). Once 
the Mott insulating plateau appears in the center of the 
system, degeneracy sets in for the lowest natural orbital 
[inset in Fig. [7|^c)], and their weight in the region with 
n = 1 vanishes [Figs.|HJe) and (f)]. 



A. Off-diagonal correlations and the lowest natural 
orbital occupation 

The peaks appearing in the momentum distribution 
function at fc = 0, and the high occupation of the low- 
est natural orbitals are consequences of the off-diagonal 
quasi-long-range correlations present in the one-particle 
density matrix (pij). We obtain that away from ni^rij = 
0, 1, Pij decays as a power law pij ~ \ixi — Xj)/a\~^^^ for 
large distances independently of the fact that the den- 
sity is changing in the system. This is shown in Fig. |U| 
where we plot pij for traps with fillings ranging over three 
decades. As an approximation, one may consider that pij 
and \{xi — Xj)/a\~^/'^ are related through a constant that 
depends only on the characteristic density p, which in 
trapped systems plays a role similar to that played by 
the density in periodic systems. 




10000 



FIG. 9: HCB one-particle density matrix vs x/a (x —\ Xi — 
Xj I) for systems with Nb = 1000, p — 2.0, rii = 0.75 (dotted 
line), Nb = 100, p = 4.47 x 10"^, Ui = 0.03 (dashed-dotted 



line), and Nb = 11, p = 2.46 x 10" 



= 2.3 X lO""* (dashed 



line). The point i was chosen in the middle of the trap and 
only j was changed. The abrupt reduction of pij occurs for 
Tij ^ 0. Thin continuous lines correspond to power laws 
sjxja. The inset shows Ao vs Nb for systems with p = 1.0 
(O)- The straight line exhibits \fNb behavior. 

A detailed study shows, however, that and \(xi — 
Xj)l a\~^l'^ arc related by a slowly varying (away from 
rii, rij = 0, 1) function of the density, which we denote as 
/2(ni,nj). We find this function to be 

h{n^.nj) ^ \n,{\ - n,)n,{\ - n,)\^l^. (24) 

In order to show it, we plot in Fig. 1101 

p^ = p,,|(x. - x,)/arV[n,(l - nO]^/^ (25) 

vs Uj for two different values of in a system where a 
Mott insulator is formed in the middle of the trap. (In 
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this case the density in the superfluid phase changes be- 
tween and 1.) The results for p'^- are then compared 
with a function yl[nj(l — nj)]^/^, where as a fit we ob- 
tain A = 0.58. The agreement between the plots shows 
that /2(?T-i, nj) in Eq. I|24(l describes very well the density 
dependence of the one-particle density matrix in harmon- 
ically trapped systems. 




FIG. 10: HCB one-particle density matrix p[ 
vs rij. p'ij is measured from points i with m = 0.5 (dashed 
line) and rii = 0.8 (dotted line) in a trap with A^i, = 951, p — 
2.61 so that a Mott insulator is formed in the middle of the 
system. The thin continuous line corresponds to 0.58[nj(l — 



The fact that the dependence of pij on the density is 
very small (a fourth root of the densities) is what allows 
one to observe the \{xi —Xj)/a\~^^^ behavior without the 
need of normalizing the one-particle density matrix. For 
very low characteristic densities {ni,nj 0) the func- 
tion /2 {rii , rij ) reduces to the exact result known for the 
continuous case /2(ni, rij) ~ {nijijY/'^ 0,|23- In Fig.lTTl 
we compare results obtained for pij and its normalized 
version, 



(26) 



The difference between the power-law decay in both cases 
is small and can be only distinguished at very large dis- 
tances. In addition, we have also plotted the results ob- 
tained for Pij normalized following KoUath et al. [35l |. 



K 
P^3 



Ptj/^rnHj. 



(27) 



As Fig. Illlshows for the HCB case, a possible power-law 
behavior in pf^ disappears much before the one observed 
in Pij without normalization. This implies that the scal- 
ing of the correlations chosen in Ref. [33 does not apply 
to very strong Hubbard repulsions. 

Considering the results for the power-law decay of the 
one-particle density matrix and for the scaling of the low- 
est natural orbital, it is possible to calculate how the 
lowest natural orbitals occupation scales in the thermo- 
dynamic limit (Aq = X^ji 4>^Pij4>^j) defined by a constant 
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FIG. 11: One-particle density matrix normalized in different 
ways (for the same system) vs a;/a (a:: =| Xi — Xj |): pij (dashed 
line), pfj as in Eq. ^'2(i^ (dashed-dotted line), and p^ as in Eq. 
I27II (dotted line). Thin continuous lines exhibit \xi — Xj\~^^'^ 
behavior. The plots were made for trap with A'';, — 951 and 
p = 2.61 so that a Mott insulator is formed in the middle of 
the system. In addition, i was fixed {m = 0.34), and j was 
changed with rij > Ui. 



characteristic density. Replacing the sums in Aq by inte- 
grals (L > a). 



Ao ~ 1/a / dx dy 

-L J-L 



cj,^ix)f2[n{x),n{y)]cj,°iy) 
|(x-y)/a|i/2 



-F{p) I I 

(28) 



where we did the change of variables x ~ XC,, y ~ Y(, 
and (/)° = R~^^-^ip^. The integral over X, Y depends only 
on the characteristic density. This is so because we have 
shown before that the scaled natural orbitals as a func- 
tion of x/(^ do not change when the characteristic density 
is kept constant. In addition, also the density profiles (as 
a function of x/Q do not change when the characteris- 
tic density is kept constant. This feature, which is valid 
for HCB's and spinless fermions, was proven in Ref. [33l |. 
Hence B2{p) and C2{p) depend only on p. The previous 
results show that with the properly defined characteris- 
tic density, the occupation of the lowest natural orbital 
scales like in the homogeneous case, proportionally to 
^/Nb- Equation H28() is valid to a good approximation for 
finite systems, as can be seen in the inset in Fig.l^where 
we plot Ao vs Nf, for different systems with characteristic 
density jO = 1, the straight line displays a \/Ni, behavior. 

A global picture of the occupation of the two low- 
est natural orbitals is shown in Fig. E| There we plot 
A' = A/ \/C,/a vs p, and compare systems with different 
curvatures of the confining potential. The Mott insula- 
tor forms in the middle of the trap when the degeneracy 
appears in the natural orbitals. The comparison between 
the plots shows that already for these system sizes (larger 
than 100 lattice sites) the finite-size corrections arc small. 
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FIG. 12: Normalized occupation of the two lowest natural 
orbitals (see text) vs p. The systems analyzed have A'' = 
1000, V2a^ = 3 X 10~^t, and occupations up to 600 particles 
(thick continuous line), and A'^ = 300, V^a^ = 3.3 x 10~*t, 
and occupations up to 270 particles (dashed line). The thin 
continuous line shows a fit of Eq. 13211 . with 62 = 1.44 and 
C2 = 0.26, to our numerical results for A'' = 1000. The inset 
displays SX'q vs N~^^'^ (see text) for traps with p = 1 (0)i 
the straight line shows the result of our fit to Eq. 1291 . 



They were found to be determined by the expression 

\o/VU~a^C2{p)-D2{p)/^a, (29) 

with D2 {p) being a function of the characteristic density. 
Equation H29() has the same form as Eq. (|18|l for the ho- 
mogeneous case, when N is substituted by C,/a. In the 
inset in Fig. [El we plot 5\'^ = C2{p) - Af, vs N-^/"^ for 
systems with p = 1. The straight line displays the result 
of our fit for C2{p = 1) and D2{p =1), and confirms the 
validity of Eq. (|29|l . In the case without the lattice [2ll |. 
the same functional form was obtained for the finite-size 
corrections in terms of the number of HCB's in harmonic 
traps. 



B. Formation of the Mott insulator 

So far we have analyzed generic features of harmon- 
ically trapped systems in a lattice. In what follows we 
study in more detail the formation of the Mott-insulating 
state in the middle of the trap. In Fig.|51we have shown 
the behavior of the one-particle density matrix in sys- 
tems with fillings below the one at which the Mott in- 
sulator appears in the middle of the trap, and once the 
Mott insulator has formed, splitting the trap in two parts. 
In all these cases power-law decays pij ^ \xi — Xj\~^^'^ 
are observed. We find that when the local density ap- 
proaches to one in the middle of the trap, the power law 
Pij ^ ~ observed in one half of the system 

starts to disappear on entering in the other half. This 
effect can be seen in Fig.iniwhere we plot [Eq. ^] 
vs \xi — Xj\ (for fixed i in one side of the trap) when the 
density in the middle of the system n 1. The density 



profiles corresponding to these density matrices are pre- 
sented in the inset. The fast decay of the one-particle 
density matrix when crossing the middle of the trap (for 
?i — > 1 in the center of the system) leads to the forma- 
tion of two independent and identical quasicondensates 
in each side of the trap. This is reflected in Fig.ll2lsince 
the occupation of the two lowest natural orbitals starts 
to be similar on approaching the formation of the Mott 
insulating state. 




1 10 100 1000 

\Xi —Xj \/a 



FIG. 13: HCB one-particle density matrix pf^ [Eq. 1261 1 on 
approaching the Mott insulating state. In the plots the trap 
curvature weis kept constant V2a? = 7.5 x 10~''^t, and the fill- 
ing was increased: Ni, = 500 (thick continuous line), A'^;, = 911 
(dashed line), N^y — 931 (dashed-doted line), and Nt — 951 
(dotted line). The corresponding density profiles are shown in 
the inset. The thin continuous line following corresponds 
to a power law 0.58([a;i — Xj\/a)^^^'^. pf^ was measured in the 
four cases displayed fixing i so that rii ~ 0.2. 

In Fig. El we have also plotted a case where the den- 
sity in the middle of the system is n ^ 0.58 (thick con- 
tinuous line). In this case one can observe very small 
deviations from the \xi — Xj]"^^^ behavior in one half of 
the system on entering in the other half. This shows that 
although the destruction of quasi-long-range correlations 
between the two halves of the trap occurs only when the 
density in the middle of the system approaches Ui = I, 
deviations from the power law [a;^ — Xj\~^^'^ are observed 
before. They just increase with increasing the density 
in the middle of the trap. Only in the case of very low 
characteristic densities is the power law \xi — Xj\^^^^ ob- 
served in the whole system as shown in Fig. 1141 A com- 
parison between Figs. 1131 and 1141 reveals an interesting 
feature, which is also present in the homogeneous case 
(Fig. 13) . The interface between the short-distance and 
large-distance behavior of the one-particle density ma- 
trix moves toward shorter distances with increasing the 
density in the system. 

The destruction of the power-law decay of pij observed 
in Fig. El is reflected by an increase of the full width at 
half maximum of the momentum distribution function 
(w) close before the Mott insulator appears in the mid- 
dle of the system, w is proportional to the inverse of the 
correlation length in the system. In the HCB case the 



10 




by Wessel et al. [36j for higher dimensional s yst ems. As 
it was shown by recent experimental results [T3 for the 
systems sizes available experimentally w gives a signal of 
when the Mott insulator appears for the soft-core boson 
case. Since local Mott-insulating and supcrfluid phases 
coexist, in order to unambiguously characterize these sys- 
tems local probes are needed. This has been discussed in 
Rcfs. |36|, ^3 fo'' soft-core bosons and Refs. [s^, 113 for 
fcrmions. The basic idea of local probes follows from the 
fact that the superfluid phases are compressible while the 
Mott insulating phases are incompressible (a gap opens 
in the charge excitations). However, such local probes 
have not been yet implemented experimentally. 
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\Xi —Xj \/a 

FIG. 14: HOB one-particle density matrix pf^ [Eq. 126ll l for 
very dilute systems with Nb = 100, p = 4.47 x 10"^ (dashed- 
dotted line), and A^i, = 11, p = 2.46 x lO"'^ (dashed line). 
The thin continuous line following pf^ corresponds to a power 



law 0.56(|a;i 



>l/«) 



-1/2 



Pij 



was measured in both cases 



displayed fixing i in one border of the trap. 



correlation length is zero in the Mott-insulating phase 
and proportional to the size of the region < n < 1 
for the superfluid phase. This implies that on filling the 
trap w starts to reduce (the system size increases) up to 
the point at which the system begins to split close before 
the formation of the Mott insulating state in the middle 
of the trap. This can be seen in Fig. El where we plot 
w as a function of the characteristic density for different 
system sizes. The points at which the Mott insulator ap- 
pears are signaled by an arrow (p ~ 2.6 2.7). The 

figure also shows that, as it can be intuitively expected, 
the increase of if becomes smaller on increasing the sys- 
tem size since the superfiuid phase (the one producing 
the peak) becomes larger. The proposal of measuring w 
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FIG. 15: Full width at half maximum of the momentum 
distribution function {w) as a function of the characteristic 
density for three different traps with ¥20^ ~ 3.0 x 10~^t 
(dotted line), Via^ = 3.3 x 10~^t (dashed-dotted line), and 
¥20^ = 3.0 x lO'^t (dashed line). The arrows signal the point 
at which the Mott insulator sets in the middle of the trap. In 
the plot w was normalized by the lattice constant a. 

in order to (approximately) detect the parameter region 
where the local Mott insulator appears in the trap was 
done by KoUath et al. [s^ for ID systems and discussed 



C. Low-density limit in the lattice 

Finally, we study the corrections introduced by the lat- 
tice to the occupation of the lowest natural orbitals in the 
confined continuous case. Due to the combination lattice- 
harmonic potential the ratio between the level spacing 
and the bandwidth in the low-energy region is propor- 
tional to (C/a)~^. As in the periodic case, we find that 
the occupation of the lowest natural orbitals approach 
their value in the continuous system A^ (Nh) in the same 
way. 



A„ [Nb, CI a) = A^ (iVb) - 



C/a ■ 



(30) 



where A,, (iV^, ^/a) is the occupation in a lattice for a 
given Nb and C, and (A^t) is a function of the number 
of particles for each 77. 

We calculate A^ (Nt) and E^-^ (A^t) for the three lowest 
natural orbitals analyzing systems with fillings up to 401 
particles and sizes up to 2000 lattice sites. In Fig. 1161 we 
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FIG. 16: Extrapolated values of the first (0)i second (A), 
and third (O) natural orbital occupations in the continuous 
system (Aq). The lines following the data for the first two 
natural orbitals correspond to the results obtained in Ref. l2ll| . 
The inset shows A A (see text) as a function of (<^/a)~^ also 
for the first (0)i second (A), and third (O) natural orbitals 
in systems with 301 HCB's and sizes between 1000 and 2000 
lattice sites, the straight lines are the result of our fits. 
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compare our results for (iV^) with the ones presented 
in Ref. for the continuous case [see Eqs. (91) and (92) 
there]. The agreement is excellent like in the periodic 
case. As an inset we plot AA = A^ (Nt,) — {Nt,, (/a) 
also for the three lowest natural orbitals as a function 
of (C/a)^^ confirming Eq. H^Ofl . (The results obtained 
for AA for the second and third natural orbitals arc very 
similar, so in the inset they arc one on top of the other.) 
We find that Eq. describes our numerical results up 
to characteristic densities very close to the the one at 
which the Mott insulator appears in the middle of the 
trap. Hence the dependence of the lowest natural orbital 
occupation on the density in Fig. 1121 can be understood 
as follows. 

Comparing Eqs. (pin|l and if^ . it is possible to deter- 
mine the dependence of Ag (iV&), [IJ and (Nf,) on Nb 

Ar, (Nb) = 62 ^/N'b, El [Nb) = nI'^ , (31) 

where 621 C2 arc parameters that can be determined nu- 
merically. Both relations above were confirmed by our 
numerical results. Then Eq. H30() may be rewritten as 

Ao(^6,p)/yc7^= V^(&2-C2 p), (32) 

which describes very well the behavior observed in Fig. 
E]beforc the Mott insulator appears in the middle of the 
system. This is shown in Fig.^jby a fit to our numerical 
results for N = 1000 with &2 = 1-44 and C2 = 0.26. 
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FIG. 17: Profiles of the two lowest natural orbitals for trapped 
systems with 700 lattice sites, Vs,a^ = 2.0 x 10~^^t, and oc- 
cupations of 201 (a),(b) and 451 (c),(d) HCB's (continuous 
line); 300 lattice sites, Vga* — 1.0 x 10~^^t, and occupations 
of 69 (a),(b) and 155 (c),(d) HCB's (dashed line). The density 
profiles have been included as thin dotted lines. 



A. Off-diagonal correlations and the lowest natural 
orbital occupation 



V. HARD-CORE BOSONS TRAPPED IN 
OTHER CONFINING POTENTIALS 

The results obtained in the previous section for HCB's 
confined in harmonic traps can be generalized for arbi- 
trary powers of the confining potential when the appro- 
priate length scale, 

C - (K.A)"'/" , (33) 

is considered, [a is the power of the confining potential 
in Eq. Q]- In this section, we perform systematically 
the generalization for a power a = 8 of the confining po- 
tential in order to show that the statement above holds. 

It has been already shown in Ref. [s^l that density 
profiles as a function of the scaled positions in the trap 
(x/C) do not change when parameters in the system are 
changed keeping the characteristic density constant. This 
holds for any power of the confining potential. We start 
this section showing that the scaling of the lowest nat- 
ural orbital, defined by [Eq. H23|) ] for the harmonic case 
(a = 2), also holds for a = 8. In Fig. El we compare the 
scaled two lowest natural orbitals for systems with dif- 
ferent filling but the same characteristic density, (p was 
defined hke in Eq. ^ with C = {Vs,/t)~'^'^ , as corre- 
sponds to a = 8. (Density profiles were included as thin 
dotted lines.) The plots show that indeed the scaling 
defined holds. 



The increase of the power a of the confining potential 
leads to the formation of more homogeneous density pro- 
files, as shown in Fig.^l The fast changes of the density 
are restricted to smaller regions in the borders of the sys- 
tem. In Fig. El it can be seen that also for these systems 
off-diagonal quasi-long-range correlations are present in 
the one-particle density matrix. Like for the harmonic 
case, they decay as power laws pij ^ \{xi — Xj)/a\~'^/'^ 
for large distances and away from ni,nj = 0, 1. In the 
figure, we show results for a very dilute case, and for a 
trap with a density Ui ~ 0.5 in the middle of the sys- 
tem. A detailed study shows that like in the harmonic 
case, pij and \{xi — Xj)/a\~^^'^ are related by a slowly 
varying (away from n^, = 0, 1) function of the density. 
We find this function, which we denote as fs(ni,nj), to 
depart at high densities from the /2("i, nj) [Eq. 126() ] be- 
havior found in the harmonic case, i.e., the functional 
dependence of pij on the density is not universal. This 
can be seen in Fig. El where we plot p [Eq. vs 
Uj for two different values of rii, and compare the re- 
sults with 0.58 [nj(l — rij)]^^^. The differences between 
the plots start to be evident for large densities. 

Considering the results above for the power-law decay 
of the one-particle density matrix and for the scaling of 
the lowest natural orbitals, it is possible to generalize 
for arbitrary powers of the confining potential the re- 
sult obtained in the previous section for the occupation 
of the lowest natural orbital in the harmonic case [Eq. 
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FIG. 18: HCB one-particle density matrix vs x/a (x —\ Xi — 
Xj I) for systems with A'';, = 400, p = 0.92, rn = 0.46 (dashed 
Une), and for a very dilute system with A^;, — 11, p = 7.56 x 
10"", rii = 1.6 X 10"^ (dotted line). The point i was chosen in 
the middle of the trap and only j was changed. The abrupt 
reduction of pij occurs for rij —> 0. Thin continuous lines 
correspond to power laws ^ xja. The inset shows Ao vs A'';, 
for systems with p = 1 (0)i the straight line exhibits 
behavior. 
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FIG. 19: HCB one-particle density matrix p\j [Eq. (I25ll 1 vs 
rij. p'ij is measured from points i with Ui = 0.51 (dashed 
line) and rij — 0.34 (dotted line) in a trap with a — 8, Nt — 
999, and p — 2.05 so that a Mott insulator is formed in the 
middle of the system. The thin continuous line corresponds 
to 0.58[nj(l - Wj)]^/", like in Fig. [Till 

(|28|l ]. The differences between fs{ni,nj), f2{ni,nj), and 
eventually an arbitrary faini, rij) do not change the final 
result in Eq. (|28|l since density profiles as a function of 
the scaled positions x/C are unchanged when the char- 
acteristic density is kept constant. Then, independently 
of the power a of the trapping potential in Eq. , the 
occupation of the lowest natural orbital is determined in 
the thermodynamic limit as 

= a (p) VN'b = D^ ip) v/c7^ , (34) 

where Ca{p) and Da{p) are functions of p for a given 
power of the confining potential a. The thermody- 
namic limit is defined keeping the characteristic density 
p = Nt/C constant with C given by Eq. H33|) . As for the 
periodic and harmonic cases, the scaling in the thermo- 
dynamic limit is valid up to a good approximation for 



finite-size systems, as is shown in the inset in Fig. 1181 for 
a = 8 and A^b > 100. 
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FIG. 20: Normalized occupation of the two lowest natural 
orbitals (see text) vs p. The systems analyzed have A'' = 
700, Vsa* = 2 X lO^^^'f, and occupations up to 501 particles 
(thick continuous line), and N = 300, Vsa* = 1.0 x 10"^^t, 
and occupations up to 181 particles (dashed line). The thin 
continuous line shows a fit of Eq. I|37|l . with fcg = 1.36 and 
cg = 0.29, to our numerical results for N — 700 (7 = 8/5 for 
a = 8). The inset displays SXq vs {C/ci)~^^^ (see text) for 
p = 1.0 (O)) the continuous line shows the result of our fits. 

A global picture of the occupation of the two lowest 
natural orbitals in the case a = 8 is shown in Fig. I20f a). 
There we plot A' = A/ C,/a vs p, like in the harmonic 
case (Fig. I12|l . and compare systems with different cur- 
vatures of the confining potential. The Mott insulator 
appears in the middle of the trap when the degeneracy of 
the natural orbitals sets in. The comparison between the 
plots shows that also for a = 8 and the system sizes cho- 
sen, finite-size corrections are small. In the following we 
check whether they have the same functional form as for 
the homogeneous and harmonically trapped cases. The 
inset in Fig.l20ls how s that this is indeed the c ase. We as- 
sume that Ao / \/C/a = Cs (p) — Ds (p) / \/Cja, like for the 
harmonic trap [Eq. H29|) ]. and calculate Cs{p) and D^{p). 
In the inset in Fig. 1201 we plot 5\'q = C%{p) — Xq/ \/C,/a 
vs {C/a)^^/'^ for p = 1, the linear behavior observed con- 
firms our assumption. 

B. Low-density limit in the lattice 

In the following, we analyze the corrections introduced 
by the lattice to the occupation of the lowest natural 
orbitals in the confined continuous case. For a trap 
with power a = 8 of the confining potential, it can be 
proven numerically (Fig. I21() that in the low-energy re- 
gion the ratio between the level spacing and the band- 
width reduces proportionally to (Vgo^/t)"'^/^, i.e., propor- 
tionally to (C/a)~®/^. As for the periodic and harmon- 
ically trapped systems, we find that the lowest natural 
orbital occupations approach their value in the continu- 
ous system A,, {Nb) following the same behavior, 

A,, (Ab, C/a) = A, {Nk) - El {N^) {(:/a)-^/\ (35) 
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FIG. 21: Level spacing (in units of t) vs Vga*/f. (O) between 
first and second level, (A) between second and third level, 
and (O) between third and fourth level. The lines following 
the data correspond to power laws (Vga^/t)"'^. 



We calculate (Ni,) and Ef-^ (iV(,) for the three lowest 
natural orbitals analyzing systems with fillings up to 401 
particles and sizes up to 2000 lattice sites. The results 
obtained arc shown in Fiss. \'2'2i a.) and l22r b'). respectively. 
In Fig. l22l we compare our results for A,, {Nb) and E^^ (Nb) 
with the power laws expected from comparing Eq. H35() 



and Eq. (PH), i.e., 

A^ (Nb) = fog ^/N'b. Es (Nb) = cg (36) 

The plots show that for large occupations, basically 
larger that 100 HCB's, the expected power laws are ob- 
served. In addition, we plot as an inset in Fig. I22f b) 
AA = A^ (Nb) — Xjj {Nb,(/o-) also for the three lowest 
natural orbitals function of (C/a)"*/^ These plots 
show that the scaling behavior proposed in Eq. H35|l is 
correct. We also find that like in the harmonic trap, Eq. 
(I35|l follows the numerical results up to fillings very close 
to the one at which the Mott insulator appears in the 
middle of the system. Equation (|35|1 can be rewritten for 
the lowest natural orbital as 

Ao {Nb, C/a) = (68 - C8 /5T) (37) 

with 7 = 8/5. Similar expressions were obtained for the 
periodic (changing p by p) [Eq. H21|) ] and harmonically 
trapped [Eq. H32|l ] systems, with 7 = 2 and 7 = 1, respec- 
tively. As is shown in Fig. EOI Eq. (with 7 = 8/5) 
describes very well our numerical results below p — 2 for 
the following fitting parameters: bs = 1-36 and cg = 0.29. 




FIG. 22: Values obtained for A,, (Nt) (a) and (Nb) (b) vs A^i,, for the first (Q), second (A), and third (O) natural orbitals. 
The lines following the data correspond to power laws x/iVt (a) and N^^''^'^ (b). The inset in (b) shows AA (see text) as a 
function of ((^/a)~*^^ also for the first (0)5 second (A), and third (O) natural orbitals in systems with 301 HCB's and sizes 
between 1000 and 2000 lattice sites. 



VI. CONCLUSIONS 



We have studied in detail ground-state properties of 
hard-core bosons confined on ID lattices. The results 
obtained for the off-diagonal behavior of the one-particle 
density matrix show that it decays with the same power 
law known from periodic systems p5l |. independently of 
the power of the confining potential. In contrast, we find 



that the small dependence of the one-particle density ma- 
trix on the density is not universal. In the harmonically 
trapped case we were able to find a function that fits very 
well this density dependence [Eq. H24|) ]. The results we 
have obtained for HCB's are equally valid for the ID XY 
model with a space- varying transverse field. 

We have shown how the occupation of the lowest nat- 
ural orbital is set by the large distance behavior of the 
one-particle density matrix. Even in the cases where a 
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region with rii = 1 builds up in the middle of the sys- 
tem, we find that this quantity scales proportionally to 
the square root of the number of particles (at constant 
characteristic density), and independently of the power 
of the confining potential. The functional form of the 
constant of proportionality, as a function of the charac- 
teristic density, was determined for fillings below the one 
at which the Mott insulating region appears in the den- 
sity profiles. We have also obtained the first correction 
to the square-root behavior of the lowest natural orbital 
occupation due to finite-size effects. We find that it is 
also universal independent of the power of the confin- 
ing potential. Although not discussed here, we should 
mention that the natural orbital occupations A,, also dis- 
play another universal behavior for large values of 77 at 
low densities ~ V~^i ^ behavior that is shared by 
the momentum distribution function for large momenta 



Finally, we studied systematically the low-density limit 
in a lattice. We have shown how the results for the oc- 
cupation of the lowest natural orbitals in continuous sys- 
tems can be obtained as an extrapolation in a lattice. 
The only input knowledge ones needs to do this extrapo- 
lation is the behavior of the ratio between the level spac- 
ing and the bandwidth for the lowest energy levels as a 
function of the strength of the confining potential. 
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